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Abstract. In this paper we study the construction of a discrete solution for a hyperbolic system of partial 
differentials of the strongly coupled type. In its construction, the discrete separation of matricial variable 
method was foUowed. Two sepárate equations in differences were obtained: a singular matricial and the 
other one a Sturm Liouville vectorial problem, which by the superposition principie yield a stable discrete 
solution. 
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Resumen. En este trabajo se construye una solución discreta a un sistema en derivadas parciales hiper- 
bólico de tipo singular fuertemente acoplado. En su construcción se siguió el método de separación de 
variables matricial discreto, obteniéndose dos ecuaciones en diferencias por separado: una matricial sin- 
gular y la otra un problema de Sturm Liouville vectorial, las cuales mediante el principio de superposición, 
producen una solución discreta estable del sistema. 

Palabras claves: Sistema hiperbólico, contorno discretizado, problema mixto y de Sturm Liouville, solución 
estable. 

1. INTRODUCCION 

Los sistemas hiperbólicos en derivadas parciales con coeficientes matriciales aparecen a menudo en 
la modelación de situaciones tales como el estudio de los de procesos de calefacción en microondas (6l 
[T3] , óptica |5 1 , cardiología 1 17| , fenómenos sísmicos en medios elásticos 1 15 | , entre otros. En este artículo 
consideraremos sistemas hiperbólicos en derivadas parciales con coeficientes matriciales constantes, 
definidos en el rectángulo [0, 1] x [O, T], los cuales caracterizan la ecuación de onda homogénea: 

(1) Euttíx,t)-Aujcxíx,t)=0, XG [0,1], 0< f < r, 

BiU{l,t) + B2Ux{l,t) =0< t<T, 



,o^ m(x,0) = /(x), XG [0,1], 

^ ^ Mf(X,0)=g(x), XG[0,1], 

donde U = [Ui,...,Umf, Utt = iUitt,...,Umtt)^ , Uxx = iUixx,-,Umxx)^; f = [fi, fmf , g = (gi, gm) ^ £ 

C", A, Al, A2, Bi, B2, E eC en su estudio es llamada condición de contorno fuertemente aco- 
plada y (|3) el problema mbcto asociado a ÍD - ©. El sistema ÍD - lE) es también llamado singular si E 
es singular de lo contraio es conciderado no singular. Una gama de casos no singulares de ÍD - © han 
sido ampliamente tratados (véase [lll, |3|, [9|, |1T1, f4l). Si bien la literatura sobre su caso singular es 
muy poca |10| estudia una mejora a [16] la cual soluciona un caso particular de {D - (|3j. Este artícu- 
lo se propone formalizar y generalizar los resultados de [lOj analizando la flexibilidad de las hipótesis 
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que se asumieron para la obtención de una solución al sistema de ecuaciones mencionado. El orden a 
desarrollar en el presente artículo es el siguiente: 

-Sección 2: Se presentan definiciones y teoremas generales de las teoría de sistemas matriciales sin- 
gulares y de las ecuaciones en diferencias finitas, en particular demostraremos 3 proposiciones indis- 
pensables para la solución de la versión discreta de O - © . 

-Sección 3: Luego de discretizar el sistemaO - ll3|, aplicando el método de separación de variables 
matricial discreto, se construye una solución a la ecuación de onda considerando su condición de con- 
torno y se enuncia el teorema que resume el resultado obtenido. 

-Sección 4: Se estudia el problema mixto y la estabilidad de la solución general. Con base en el teorema 
anterior enunciamos el segundo que garantiza la existencia de una solución discreta y estable de (D - . 

-Finalmente complementaremos lo obtenido con un ejemplo ilustrativo y enunciaremos algunas ob- 
servaciones y conclusión finales. 

2. PRELIMINARES 

En esta sección damos algunos resultados y notaciones que serán de utilidad para las siguientes 

Definición 1. SeaAe c'"'"", una matriz, la cual denotamos por A'^, es llamada inversa generalizada de 
A si verifica la condición 

AA^A = A 

Teorema 2. {Teorema de Mitra) Sea A e c'"'"" y A'^ una inversa generalizada para A. Sea b e C", 
entonces 

Ax = b tiene solución o AA^b = b 
y cualquiera de sus soluciones tiene la forma 

(4) x^A^b + [l-A^A)z, zgC". 

Demostración: Véase 112] . 

Definición 3. Sea W un subespacio de y A una matriz en c'"'"", se dice que W es invariante por la 
matriz A, si verifica que AW c W {i.e ^xeW ^ AxeW). 

Teorema 4. Sean A, B, B'^ e c'"'"", entonces B A [l - B'^ B] = O siiKer {B) es un subespacio invariante de 
A. 

Demostración: Véase 112] . 

Definición 5. (La Inversa Drazin) Sea A £ c'"'"" tal que ind{A) - k, la inversa Drazin de A (A^) es la 
única matriz que cumple 

1) A^AA^ = A^ 

2) A'^A = AA^ 

Teorema 6. Este teorema consta de 3 numerales: 
1. Si A tiene la descomposición canónica de Jordán 

I C Q 



(5) A=T 



O N 



donde C e CP^'p yTe c'"'"" son matrices invertibles, N e C^""^ es una matriz nilpotente con índice k con 
[p, q) en cumpliendo p + q-n, entonces la inversa Drazin de A, aP , se obtiene de la siguiente forma 



(6) A^ = T 



c-i O \ . 



O O 



T~ 
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2. Para cierto polinomio Q (x) , La inversa Drazin de A puede ser escrita como 

= Q [A) 



' c-^ O ^ 



3. Si en la descomposición canónica de A en N -O, entonces 

:-i o 
o o 

Demostración: Véase H] y [12] . 

Definición 7. Sea A e C^""^, / : C ^ C, decimos que f pertece al espacio 9 {A) si existe un entorno V del 
espectro de A, a {A) sobre el cual fes analítica. 

Teoremas. [Aplicación espectral) Seafe entonces 

fia{A))=a{f{A)], 

donde f {a {A)) = {/ (A) : A g o- (A) } . 
Demostración: Véase (3 . 

m 

Definición 9. Sea x e C"^ . definimos la norma vectorial \\\\i por; ||x||i = ||(xi,...,Xto)IIi = Z \Xi\,ylanorma 

i=l 

matricial por ella inducida como 

(7) llalli =sup{||Ax||i:xeC'", ||x||i = l}, paraAeC""""" 

si las columnas de A son ai, az,—, a„ para las cuales definimos sus componentes aij, i - l,..., n, para 
i - 1, n respectivamente, entonces l|3 cohincide con 



llalli =máx|^^ \aij\:j = 



Teorema 10. Sea A e c'"'"", / e 3 (A) , la función matricial f (A) se define por la matriz. 

f{A) = P{A) 

donde P(x) es un polinomio conocido de grado menor que el del polinomio minimal de A. 
Demostración: Véase |7|. 

Teorema 11. Sea Ae C^""" talque a - max {\z\; z e a (A)} , para e > O existe una norma matricial en 
C'"^'" talque \\A\\<a + e. 

Dada la variedad de tipos de estabilidad manejadas en la literatura de los esquemas de diferencias fi- 
nitos numéricos asociados a ecuaciones en derivadas parciales, nosotros hemos optado por la siguiente 
definición: 

Definición 12. \l^{Estabilidad). Consideremos fijo el par ordenado {h, T) eW^ y. . Si se tiene una 
malla de puntos en el plano (x, t) con paso espacial h y paso temporal k > O talque M es un número 
natural. Si U es una función definida en dicha malla y U[ih, jk) representa su valor en el nodo [i, j) , se 
dice que U es estable, cuando U[ih, jk) permanece acotada independiente de los valores dekyM tal que 
kM= T 

Definición 13. (Problema de Sturm Liouville Discreto) El problema de valor en la frontera formado por 
la ecuación en diferencias 

(8) ^[p{i-\))^u{i-\) + q{i)u{i)+Xr{i)u{i)=Q, 0<i<N 
con condiciones de frontera 

(9) u(.0) = aua), um^ PuiN-l), 
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en las que las funciones son escalares se llama problema de Sturm-Liouville discreto (PSD). 
Donde A es el operador de diferencia definido por 

(¿ -D^p (/) -p{i-l), A^pm = A [Ap [i]] . 

en la ecuación en diferencias Á.esun parámetro, a y ¡5 son constantes conocidas, la función p (/) esta 
definida enO< i < N-l, p (/) , qii) y r {i) enO< i < N; además pii)>0 yr {i) > 0. 
Al efectuar la operación A en la ecuación en diferencia^, encontramos para 0< i < N 

p (/) u (¿ + l)-[p ü) + p{i- 1)) u (/) + [q (i) + Ar (/)) u (¿) +p{i-l)u (¿ - 1) = O 

denotando sii) = pii) + pii -l) - qii) ,0 < i < N, se tiene 

(10) -p{i- l)uii - 1) + 5(¿) u{i) -p{i)u{i + 1)- Ar (/) u{i) , O < i < N 
en esta última ecuación para los dos casos k=\,N - \ tenemos las ecuaciones 

-pmu{Q) + s{\)u{\)-p{\)u{2) = Ar(l)M(l), 
-p{N-2)u{N-2) + s[N-l)u{N-l)-p[N-l)u[N) = Ar (Ai- 1) M (Ai- 1) , 
las cuales teniendo en cuenta @ se reescriben como 

(11) smuW-pWuiZ) = Ar(l)u(l), 

(12) -piN-2)uiN-2) + siN-l)uiN-l) = XriN-l)uiN-l), 
dondes{l) = sa)-ap{.0)ys{N-l) = s{N-l)-Pp{N-l). 

Con las ecuaciones {TTJ para i = O, 03 para i = N- 1 y 03 para 2<i< N-2 se producen N- 1 ecuaciones 
equivalentes a @ - @ , te cuales describen un problema generalizado de valores propios matricial con la 
forma siguiente 

uW 
u(2) 



u{N-2) 
[ u{N-l) 

uH) 
u{2) 

u{N-2) 
u{N-\) 

Ya que p {i) > O, O < i < N - l; r (i) > O, O < i < N de acuerdo a la teoría espectral para matrices simétricas 
nosotros asumimos los siguientes teoremas. 

Teorema 14. El PSD tiene exactamente N-l eigenvalores A/, O < Z < N; a cada uno de los cuales les 
corresponde las eigenfunciones vi{i), 0< i < N definidas para 0< i < N. 
Demostración: Véase UJ ó 114] . 

TeoremalS. Sean Xi y viii) losN-l eigenvalores y egenfuciones del PSD, entonces el conjunto 

{viii):l = l,...,N-l} 
es ortogonal con respecto al producto interno 

N-l 

Y, r (í) v^ (¿) Vy (¿) =0, Vju 7^ V. 



(13) 



5(1) 

o 

O 



-pW 

5(2) 
-p(2) 



O 



r(l) 
O 

O 
O 



o 

r(2) 
O 



siN-3) 
-p{N-2) 



-p{N- 
s{N-\ 
O 



O 

r{N-2) O 
O r(Ar-l) 
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Demostración: Véase fTl . 

Teorema 16. Sea {vi{i) :0 < l < N} el conjunto de eigenfunciones en el teorema anterior y u (¿) cualquier 
función definida para O < i < N , entonces u (/) puede expresarse como una combinación lineal de dicho 
conjunto, es decir: 

N-l 

u Ü) = ^ Q i;; (í) , 0<i<N 



donde cada ci puede ser determinada por 

H^-^rii) vi{i)u{i) 

ci = 



0<1<N 



Lfj.'rii) v]{i) 
Demostración: Véase H] . 

Las siguientes 3 proposiciones las enunciamos con una notación pertinente para la demostración de nues- 
tro teorema fundamental Í22^ . 



Proposición 17. -.Sean ay escalares conocidos, N eNyX un eigenvalor del PSD 

(14) h[i + l)-{2-X)h{i) + h{i-l) = 0,0<i<N, 

(15) h{.0) + aNiha)-hm = O, 

(16) hm + PN{hm-h{N-l)) = O, 

con eigenfunción v{i). Entonces v {i) es solución del PSD - (O si y sólo si 

H{i) = v{i)R 

con R vector arbitrario en C", es una solución del problema de Sturm Liouville discreto vectorial (PSDV) 

(17) H(i + l)-(2-A)H(z) + H(z-l) = 0,0<i<N, 

(18) H{0) + aN{Ha)- Hm 

(19) HiN)+l3N{H{N]-HiN-l)) 

en el que H{i) es una función vectorial discreta en C™. 
Demostración: 

Supongamos que H (i) satisface el PSDV Ü^-Ü^ dondeH(i) - ( h^ (i) h2{i) ... hmíi) ] siendo hic {i) 
funciones escalares discretas definidas para O < i < N. Luego en virtud de la igualdad entre vectores tene- 
mos que cada una de sus componentes h]c (/) A; = 1, m satisfacen el PSD lO - fZ^ . Entonces al hallar 
una solución para una componente, ella lo será también para las demás. Sin pérdida de generalidad de- 
notemos hic {i) = h {i) . Si en lO definimos: 

q{i)=0, r{i) = 1, p{i) = 1, 

de 03 - 03 se obtiene 

■LiML -10 O 
-12-1 : 



O, 
O, 



(20) 



O 



: 2 -1 

O -1 ^ 

Por el teorema lO , l[20l tiene los N-l eigenvalores 



■ h{l) ■ 




' h{\) ' 


h[2) 




h{2) 




= X 




h{N-2) 




h{N-2) 


h{N-\) 




h{N-D 



(21) 



{Xi 



N-l 
'1=1 
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con eigenfunciones asociadas 

(22) {i^z(¿)}{l"/ , 0<1<N, 0<i<N, 

vi^ivdl), vim,..vdN-l))^, 0<1<N. 

Luego, las soluciones del PSD lO - ([7^ tienen la forma 

vi{i) - Vi (i) D2, con D2 arbitrario enC e i = l, ...N - 1. 
Como lo anterior se verifica para cada componente de H{i) concluimos que 

(23) H{i) = {hi (/) , hm (¿)) = Vi (í) R, R arbitrario en C™ 
es solución del PSDVÜT^ - 03 . 

Aplicando el teorema {73 es claro afirmar que el PSDVÍT?^ - lO tiene exactamente N -l soluciones de 
la forma C23II donde Vi {i) son las N-l eigenfunciones de PSD - iO determinadas en Í2^ . 

Proposición 18. : SeaH{i) función vectorial discreta en C", que verifica 

(24) H{i + \)-{2-X)H{i) + H{i-\) = 0,0<i<N, 

(25) AiH{0) + A2NiHa)-Hm = O, 

(26) BiHiN)+B2NiHiN)-HiN-l)) = O, 

con Ai,A2,Bi,B2eC"""^y 

tal que rang[G[a, P)) < myRe KerG[a, f)). Entonces para cada Xi,0 < l < N eigenvalor del PSD 03 - 

Hiü) = vi{i)R, 0<i<N 

son N-l soluciones de |[24t - (|25t . 
Demostración: Sea 

Hiü) = viii)R0<i<N,ReKerG{p,l5] 
con H{i) así definida, - (|25t es equivalente a 

(28) ivii + l)-{2-X)vii) + vii-l))R = 0,0<i<N, 

(29) iAiv{0) + A2N [vil) -vm)R = O, 

(30) (Bii;(Ar) + B2A^(í^(Ar)-í^ (Ai -!)))« = O, 

como V (z) verifica C3 - (O , se tiene que (2^ - , es equivalente a 

{v{i + Y)-{2-X)v{i) + v{i-Y))R = 0,0<i<N, 
(ayli -A2) (1/(1)- 1/(0)) = O, 
[PBi-B2]iviN)-viN-l))R = 0. 

Como además ReKerG[a,l3], se concluye que Hi{i) - vi (z) R, satisface (2^ - lH^ . 

El teorema y proposición siguientes los enunciamos con una notación pertinente para la demostración de 

nuestro teorema fundamental. 

Teorema 19. Sean I, Ae c"""" donde I es la identidad, yG[j) función vectorial definida para i = O, ...N; 
supongamos además p escalar y consideremos la siguiente ecuación en diferencias matricial 



G (7 + 1) - (27 + pA) G [i) + G (;• - 1) = O, ; > 0. G (;) e C 



mxm 
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Entonces su solución general es 

G[j) = [P+iA)]jPi + [P-iA)]jQi, j>0, h,l2 gC", 

donde P+ (x) y P_ (x) son los polinomios de grado más pequeño que el grado del polinomio minimal de A 
tal que 

P+ {A) = 7+|a+^(/+|^)'-7,P_ {A) = I+^A-yJ[l+^Af-I. 
Demostración: Véase 13] , o también HU . 

Proposición 20. Con las mismas hipótesis del teorema anterior, consideremos la siguiente ecuación en 
diferencias matricial 

(31) EG[i + \)-{2E + pA)G[i) + EG[i-\]^i), ; > 0. 

donde i? e C'"'"" es singular, supongamos además que existe jeC tal que {jE + A) es no singular, defini- 
namosE = [jE + A)~^E, A = [XE + A)~^ A. Entonces su solución general viene dada por 

G[j) = Z^^ ÉÉ^h + Z^^ÉÉ^h, con h, h arbitrarios en C"" . 
donde Zq = [P+ [É^Á]] ÉÉ^, Z^ = [P- [É^ Á)] ÉÉ^ 

con 

P+ [É^Á] = 1 + ^É^Á + yJl^I+^É^Á^^ -I, P- [É^Á] = 1+ ^É^Á - 

Demostración: Dado que existe y con [yE + A) invertible, sean É = [jE + A)~^E, Á= {XE + A)~^A, luego la 
ecuación Í3]\ es equivalente a 

(32) ÉG[i + \)-{2É-pÁ)G[i] + ÉG[i-\) = Q, j>0. 
Considerando la descomposición canónica de É, obtenemos: 

en lasqueCeCP^P es matriz invertible, N e C'''' es una matriz nilpotente de índice k con p, q cumpliendo 
p + q - n. Luego haciendo 

i Hj) 



[ dU) j' 



(34) TG[j) 
Í31] es equivalente a los dos siguientes sistemas: 

(35) Ch[i + \)-[2C + p[l-YC]]h[i) + Ch[i-\]^{), 

(36) Nd[i + \)-[2N + p[l-YN]]d[i) + Nd[i-\] = {). 
Como C es invertible se obtiene que lO se verifica si y sólo si 

h [j + 1) - [21 + pC-' [I -rC]]h [j] + h[j-l] = o, 

la cual por el teoremalT^ tiene la solución general 

h [j] = [P+ iD)]j Pi + [P_ {D)]j Qi, i > O, h, h e CK 

donde D = C~^ [l - jC) y P+ (x) y P- (x) son los polinomios de grado más pequeño que el grado del poli- 
nomio minimal de D tal que 



PAD) = I+^D + J(l + ^Df-I,P.{D) = I+^D-J(l + ^Df-I 
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Como para N suficiente mente grande 2N + p[l - jN] es invertihle y k es el Índice de nilpotencia de N, 
por inducción puede demostrarse que la ecuación tiene la solución trivial d [j] =0, j > O, luego con- 
siderando 

,37, n)^-^'-^-ÍVI]^ 

y haciendo uso de las ecuaciones lO , 03 y (3B , se concluye que tiene la solución general 

(38) G[j] = Z^^ ÉÉ^h + Z^^ÉÉ^k, con h, h arbitrarios en C"". 

(39) donde Zq = [P+ [É^Á]] ÉÉ^, Z^ = [P_ [É^Á]] ÉÉ^ 

3. Solución discretizada de la ecuación de onda 
Retomemos nuestro problema central 

(40) Euttíx,t)-AUxAx,t) = O, XG[0,1], í>0, 

(41) AiUÍ0,t) + A2Uxi0,t) = 0,t>0, 

(42) Biua,t)+B2uAl,t) = 0,t>0, 

(43) u{x,0) = f{x), xe[0,l], 

(44) utix,0) = g(x), xe [0,1], 

En el que asumimos además la hipótesis: que existe y g C tal que {jE + A) es no singular, y así definimos 
É ^ijE + Ay'^E , Á^UE + A)-'^ A conldL condición: 

(45) z^O para algún zea [É^Á] c C 

En el esquema de diferencias correspondiente a Í40[ - |[44ll dividimos el dominio [0, 1] x ]0, oo[ en rectán- 
gulos de lados Ax - h, At = k, luego al introducir coordenadas de un punto típico de la malla [ih,jk) , 
representamos el valor u[ih,jk) por U[i,j). Si aproximamos las derivadas parciales de u utilizando 
diferencias avanzadas para las derivadas primeras y diferencias centrales para las derivadas segundas. 
Ellas tienen la forma 



(46) Ut{ih,jk) 

(47) Ux{ih,jk] 

(48) Uxx[ih,jk] 

(49) Utt[ih,jk] 



U[i,j + l)-U[i,j) 
k 

U[i + l,j)-U[i,j) 
h 

U[i + l,j)-2U{i,j) + U{i-l,j] 

U{i,j + l)-2U{i,j) + U{i,j-l) 
A;2 



Al sustituir |[46ll - 1[49] | en |[4Dl l - (l44l l , para N un entero positivo con h = l/NYr = k/hO<i<N;se obtiene 
la siguiente representación aproximada de (|4ÜJ - (l44l l 



(50) 



r^A{U{i + l,j)-U{i-l,j))+2{l-r^Á]U{i,j)~ 
{U{i,j + l) + U{i,j -l)) = 0, 0<i<N, j>0 
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(51) Aimo,j)+NA2[ua,j)-mo,j)] = o, j>o, 

(52) BimN,j) + NB2[mN,j)-mN-l,j)) = 0, j>0, 

(53) Uii,0)^F{i),0<i<N, 

Uíi,l)-U{i,0) ^^.^ ^ . 

(54) ; = Gii), 0<i<n. 

k 

Busquemos soluciones no triviales a la ecuación de onda discretizada llHÜll . Para esto usamos el método 
de separación de variables discreto matricial suponiendo U[i, j) de la forma 

(55) U[i,i)^G[i]H{i), H{i)e£"'. 

Para r - k/h, imponiendo que {U[i,j)} verifique ÍHÜ), resulta 

r^AG{j][H{i + l)-H{i-l,j)]+2{E-Ar^)G[j]H[i)- 
^ ^ E[G{j + l) + U{j-l)H{i)] =0,0<i<N, j >0 

Tomando peU. arbitrario y puesto que 2[E- Ar^] = - r^yl ^2 + j + [2E + pA) , |[56ll toma la forma 

^3-,^ r^AG[j)[H{i + l)-A[2 + -^]Hii) + Hii-l) - 

[EG{j + l) - {.2E + pAG{j) + EG{j - l)] H ü) = O, 0<i<N, j>0. 

La ecuación ^} se verifica si {Híí)}y {G[j]} satisfacen las ecuaciones en diferencias 

(58) H(/ + 1)-(2 + -^)h(/) + H(/-1) =0, 0<i<N, 

(59) EG[j + l)-i2E + pA)G[j) + EG[j-l] = 0, j>0. 

Si en (l58j definimos p = -r^A; para cada A; eigenvalor del PSD {TU - {TBli , de la proposición {TTl , sabe- 
mos que si Vi (¿) es el eigenvector asociado a A/. Entonces: Hi (i) definida como 

(60) Hi (i) = Vi (/) i?, i? e - {0}, 0<i<N, j> 0. 
satisface ||58). 

Por otra parte usando la propocisión GOJ, para p = -r^Xi, la ecuación ll59), tiene solución dada por 

(61) G{j) = Zq-' ££^/i + Z/££^/2, con /i, fe arbitrarios en C'". 

(62) donde Zq = [P+ [É^Á]] ÉÉ^, = [P. [É^Á]] ÉÉ^ 
así, una solución de llHÜl es de la forma 

(63) U[i, i) = zJHu (/) + Zo^H2z (¿) , O < / < AT, 7 > 0. 
donde 

(64) Hu (¿) =Piv (/) , H21 (/) = QiV (¿) con Pi, QieC^- {0} . 

3.1. Condición de contorno discretizada. Las soluciones U[i, j) de llHÜll encontradas en |[63), se escri- 
ben de acuerdo a ll64| . como sigue 

(65) Ui [i, j) = zJPi Vi (/) + Z/ Q/ Vi (í) , O < ¿ < iV, 7 > 0. 

donde vi (/) es alguna eigenfunción del PSD ÍHl - tUl y P;, Q/ e (C" - {0}) . Sea G [a, definida en ll271l , 
por la proposición {TBJ , dejando ; fijo, llüHl satisface dSD - ll52]| si Pi y Q/ en los sistemas de ecuaciones 

Z^^hAiVÍ0) + A2Niva)-vi0)))Pi = O, 

Zj iBiVÍN) + B2NiviN)-v{N-W)Pi = O, 
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y 

Z^j iAiVÍ0) + A2Niva)-vm))Qi - o, 

Z^^ iBiVÍN)+B2NiviN)-v{N-mQi = O, 

pertenecen a. Ker G [fj., p]n{,C'" - {0}). Con lo anterior [7/ (/, ;) en ll65) es solución de llHD) 
Los resultados obtenidos hasta el momento, los resumimos en el siguiente teorema. 

Teorema 21. Supongamos f{x), g(x) funciones en C"^, A, Ai, A2, Bi, B2; E e c™'"", E matriz singular 
con Y tal que (jE + A) es invertible, A, E, E^ ^ £ mxm dgjjYiidas por lO, Í33 y que verifican la condi- 
ción lO; sean además (a,/3) e que verifican lO - 03. G(a,j6) representada por ll23 que cumple 
rang [G [a, fi)) < m. Entonces una solución de ^5D^ - (5^ esta dada por: 

N-l . 

(66) U{i,j) = X [z^ÉÉ^Pi + ZlÉÉ^Qi] Vi {i) ,0<i<N, j>0 

1=1 

donde Zq y Zi están definidas por (6^ , /i, ¿2 ^ KerG[p,,p) n (C" - {0}) y vi {i) son las N-l eigenfunciones 
asociadas al PSD QH - 03 determinadas según la proposición G3 ■ 

4. Problema mixto 

Analizamos ahora las condiciones a imponer para que U[i,j) en (l66l l verifique el problema mixto 
(|53] | - (IMll • Para esto al sustituir |[66li en dichas expresiones, se deben cumplir las siguientes dos ecuacio- 
nes: 

N-l 

(67) FU) = ^ ( ÉÉ'^Pi + ÉÉ^Qi] VI (/) 



(68) G(/) = 



i=\ 

I^-i [Z^ÉÉ^Pi + Z.ÉÉ^Qi) Vi ii)-F (/) 



k 

con lo anterior, según ÍHli , cada componente de F (¿) y G (/) tienen su representación discreta respecto 
a {vi {i)}fSi^, luego las ecuaciones llBTIl y (l68) son equivalentes a 



(69) EE'^Pi + EE'^Qi 



(70) ZoEE^^Pi + ZiEE^^Qi 



que es un sistema matricial cuyas incógnitas son los vectores Pi y Q/ los cuales son soluciones de las 



ecuaciones 



^71) \Zi-Z^)EE Pi = — — 



i72) [Zi-Zq)ÉÉ^Q 



viUY 

t^Dr, _ lf.V vdi)[[Zo-l)Fii)-kGii)] 
Lfji'viiif 



Sin pérdida de generalidad, suponiendo que a (D) no contiene a O, por el teorema (|2ÜJ , se tiene que 

aiD)=a [É^A] - {0} , además 

(73) a (P+ (D)) = a (Zq) - {0} , o" (P+ (D)) = a (Zi) - {0} 



SOLUCIONES DISCRETAS PARA SISTEMAS MATRICIALES EN DERIVADAS PARCIALES. 



11 



O sea 



luego 



(7(P±(D)) = 1 1 + |d ± ^(i + ^df -l:dea [ÉÁ^] - {0}| 



o-(P-(D)-P+ (D)) = |2^(/ + |d)'-l:deo-(EÍ^)| 

y así P- (D) - P+ (D) es invertible si se verifica 

(74) prifl + ^dj^^O 

Luego tomando en |[74li p suficientemente pequeño, utilizando la inversa de P_ (D) - P+ (D) se pueden 
emplear los numerales 1 y 3 del teorema © , para hallar ((Zi - Zq) ÉÉ^) ^, la cual verifica: 

((Zi - Zo) ((Zi - Zo) = ÉÉ^ 

Así, usando el numeral 2 del teorema ll2j , se sigue que ÉÉ^ conmuta con Zq y Z^ y por lo tanto los 
sistemas JTTJ y ||72j son consistentes si cumplen la condición 

(75) ÉÉ^F{i) ^F{i), ÉÉ^GU) = G(/) 

teniendo presente lo anterior, por el teorema 0, las ecuaciones (|7D y (|72j admiten las siguientes solu- 
ciones 

^76) Pi^[[Z,-Z^]EE \ , 

L/=i viW 

^77) Qi = [[Zi-Z,)EE J . 

L¿=i i^z(í) 

Con las condiciones 

(78) {G(/),F(/), 0<i<N}czKerG{^i,l3], 

G[^,l3]É^Á[l-G[^,l3yG[^,P)]=0, 

por el teorema (m, KerG[fi,f}) es un subespacio invariante de (E^A). Por lo tanto con ||75j Pi y Q/ 
cumplen 

M^Pz = Pz, ÉÉ^Qi ^Qi, {Pi,Qi,0<i<N}^ KerG (/i, ¿6) . 

lo que garantiza la consistencia de la condición de contorno discretizada (ISTT l - Í52\ junto con el proble- 
ma mixto ll53j - (IMll • Luego, una solución de llHDJ - (I54ll viene dada por 



(79) U [i, j] = Y. [z^EÉ^Pi + ZlÉÉ^Qi] vi (/) , 



JV-l 

E 

Z=l 

con Pi y Qi determinados por (ITüll y G?} . 
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4.1. Estabilidad. Para el análisis de la estabilidad de la solución ll79), suponemos que -F (7) y G[j) son 
acotadas. Por los teoremas (dü) y ll8), de (173), el espectro de las matrices Zg y Z^ en l[79) comprende los 
valores 



dado que 



I+^d + J{l + ^df-I 



l,\/dea [ÉÁ^], considerando que p puede verificar 
p < max{\Xi\; A/, 1 < / < N- 1} x r^, 



(80) 

por teorema (HU se obtiene 

(81) ¡Zol < ||ZJ < l + 0(A;),||Zo-/ I < ||Zi-/ I < 0(A;) 

||(Zo-Zi)ÉÉ^|| < O(fc-i), A;-0. 

para alguna norma matricial || || en c'"'"". Usando (TTIT l , l|72j, (ISTT l y el acotamiento de F{i) y G(/), con- 
cluimos que 



(82) 



|P/II = 0(1), ||Q/|| = 0(1). 



Z^ 
^0 



permane- 



De lo anterior, se sigue que {U[i, j]} en |[79) permanece acotada, si los números 

cen acotados cuando ; ^ 00, A; ^ O, 1 < j < M, Mk - T. Nótese que de || Zq || < 1 + O (A;) se obtiene que 
para alguna constante positiva S, || Zq || < 1 + kS, entonces, para 1 < j <M se calcula 



(83) 



Z^ 
^0 



M ^ „MO{k) ^ „MkS - gTS 



I Zo I ' < (1 + O (A;))^ < (1 + O {k)r < e^""^"^ < e 



. Así, considerando ll82l, ll83l y L = máx{|| i;; || : 1 < Z < iV- 1} , se sigue que la 



lo mismo ocurre para 

solución definida por ÍZHJ , llZHli y llTTI es estable, es decir: 

I [7 (¿,7) I = 0(1) , A; ^ O, h = — paso temporal fijo 
1 < i < N-l j ^00, con Mk = T 
En resumen, el siguiente resultado ha sido establecido. 

Teorema 22. Con las mismas consideraciones del teorema anterior, sean G{i) yF{ i) acotadas que satisfa- 
cen ^ZB y con p verificando y (8Q ; entonces una solución estable del problema discreto (EOi - (53 
esta dada, por 



(84) 



N-l . 

U[i,j] = Y. \Z¿ÉÉ''Pi + ZlÉÉ''Qi\vi{i) 
1=1 ^ ' 



donde Pi, Qi son como lo indican las ecuaciones y ÍZ3 
4.2. Ejemplo. Consideremos las matrices en C^""^ 

e O 6 
O 7 O 





í ° 





^ 




, A = 





6 





, {e, Y, or} eC-{0} 


j 


l 





0- ; 





' (^w <^\2 b\ ^ 
azi ^22 b2 
a^i ^32 bs j 



,A2 



juaii /iai2 ci 

p.a2l /ifl22 C2 



. B2 = TfBi 
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( Mi) ] 




' gi U) ' 


FU) = 




, Gii) = 






[ J 




[ ¡ 



escogiendo A 7^ O tal que Ay + 5 7^ O, se calcula 



f 1/A O 5/Xe \ 

O r/(Ar+5) o 



( l/Ae O -6/ea ' 

O l/(Ar + 5) O 

O O 1/a j 

(O O 

, Á = 



O 5/(Ar + 5) 



-6/e \ 
O 







l 









; 


l 







1 


1 










í 









' 1 6/e ' 






í ° 










É 


D _ 


(Ar 


+ 5) /y 





,ÉÉ^ = 


1 


,É^Á = 





5/7 











l 





J 




,0 , 






l 








J 


si en ÍHli - {TBl escogemos, a 




= T], entonces 


















í ° 








^ 


T 








( 






























l/{dl + dl + 


di) 



















, di 


<¿2 ^3 





i 










, di 


d2 











donde di = ¡j.bi - c¿, / = 1, 2,3 y se comprueba además que: 

a[É^ i) = {0, 1, 3},ran(G(2,l)) <3, G(a,/3) = O, G(a,)6) G(/) = O, 
ÉÉ^F{i)^F{i), ÉÉ^G{i) = G{i), G[yi,0)É^Á[l-G[yi,0)^ G[^,p)^=0. 

Luego, si p es suficientemente pequeño, las hipótesis del teorema ES pueden ser verificadas y por lo 
tanto una solución de llHO) - es de la forma 



N-l . 

U[i,i)^ E \Z¿ÉÉ''Pi + ZlÉÉ''Qi\vi{i) 
1=1 



deonde Pi y Qi son como lo indican las ecuaciones llTBJ y ^7} y vi (z) son las eigenfunciones descritas 
en la proposición ^7} . 

4.3. Observación. 

1 . El procedimiento para el hallazgo de la solución del problema (|4DJ - imi aquí estudiado, incluye el 
caso de su versión no singular. En efecto, sin perder generalidad, supongamos que en la ecuación 

' K O 



O N 



J r-\ la matriz M + A 



T ^ Escogiendo A suficientemente grande, esto 



(So) , £■ = 7. Al considerar la descomposición canónica de A - T 

. . A/ + 7: O 

obtiene la forma a7 + = r „ t 

[ O M + N 

hará que la diagonal de los bloques A/+iCyA/+Aínose anule y así, la condición: que existe 7 e C 
tal que [jE + A) es no singular, es verificada, y nuestro procedimiento puede ser aplicado. Nótese 
además que E E = ly E A = A, por lo tanto la condición i75t puede omitirse. 
2. Si la condición ||45j del teorema |22] no se cumple, entonces de acuerdo a la proposición 123 la 
ecuación ||59j puede ser reducida a lEHJ ; como la solución de esta es la solución trivial, se concluye 
que así también lo será la solución encontrada para ll5ÜJ - l|54ll . 
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4.4. Conclusión. En este trabajo se ha demostrado la existencia de soluciones numéricas estables pa- 
ra una gama de sistemas de ecuaciones en derivadas parciales singulares de la ecuación de onda ho- 
mogénea como la descrita en iHDll - iSDi. En nuestra metodología, con respecto a los eigenvalores del 
problema de Sturm Liouville discreto asociados a lHH - ím, hemos asumido el teorema espectral para 
matrices simétricas en vez de aplicar métodos numéricos para establecer su existencia. Dicho enfoque 
nos permitió enfrentar con otra estrategia el problema de encontrar tal solución a nuestro sistema de 
interés para el cual la arbitrariedad de los parámetros a y /3, relacionados al PSD {T4J - {TBJ puede ser 
usada con el fin de que las hipótesis del teorema (|22j se verifiquen. 
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